Wetting on a spherical wall: influence of liquid-gas interfacial properties 
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We study the equilibrium of a liquid film on an attractive spherical substrate for an intermolecular 
interaction model exhibiting both fluid-fluid and fluid-wall long-range forces. We first re-examine the 
wetting properties of the model in the zero-curvature limit, i.e. for a planar wall, using an effective 
interfacial Hamiltonian approach in the framework of the well known sharp-kink approximation 
(SKA). We obtain very good agreement with a mean-field density functional theory (DFT), fully 
justifying the use of SKA in this limit. We then turn our attention to substrates of finite curvature 
and appropriately modify the so called soft-interface approximation (SIA) originally formulated 
by Napiorkowski and Dietrich [Phys. Rev. B 34, 6469, (1986)] for critical wetting on a planar 
wall. A detailed asymptotic analysis of SIA confirms the SKA functional form for the film growth. 
However, in this functional form SKA approximates surface tension with that of a sharp interface. 
This overestimates the liquid-gas surface tension and thus SKA is only qualitative rather than 
quantitative. On the other hand, by relaxing the assumption of a sharp interface, with e.g. even 
a simple "smoothing" of the density profile there, improves the predictive capability of the theory 
markedly, making it quantitative and showing that the liquid-gas surface tension plays a crucial role 
when describing wetting on a curved substrate. In addition, we show that in contrast to SKA, SIA 
predicts the expected mean-field critical exponent of the liquid-gas surface tension. 

PACS numbers: 05.20.Jj, 71.15.Mb, 68. 08. Be, 05.70.Np 



I. INTRODUCTION 

The behavior of fluids in confined geometries, in par- 
ticular in the vicinity of solid substrates, and associated 
wetting phenomena are of paramount significance in nu- 
merous technological applications and natural phenom- 
ena. Wetting is also central in several fields, from en- 
gineering and materials science to chemistry and biol- 
ogy. As a consequence, it has received considerable at- 
tention, both experimentally and theoretically for several 
decades. Detailed and comprehensive reviews are given 
in Refs. 

Once a substrate (e.g. a solid wall) is brought into 
contact with a gas, the substrate-fluid attractive forces 
cause adsorption of some of the fluid molecules on the 
substrate surface, such that at least a microscopically 
thin liquid film forms on the surface. The interplay be- 
tween the fluid-fiuid interaction (cohesion) and the fiuid- 
wall interaction (adhesion) then determines a particular 
wetting state of the system. This state can be quantified 
by the contact angle at which the liquid-gas interface 
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meets the substrate. If the contact angle is non-zero, 
i.e. a spherical cap of the liquid is formed on the sub- 
strate, the surface is called partially wet. In the regime 
of partial wetting, the cap is surrounded by a thin layer 
of adsorbed fluid which is of molecular dimension. Upon 
approaching the critical temperature, the contact angle 
continuously decreases and eventually vanishes. Beyond 
this wetting temperature one speaks of complete wetting 
and the film thickness becomes of macroscopic dimension. 
The transition between the two regimes can be qualita- 
tively distinguished by the rate of the disappearance of 
the contact angle, which is discontinuous in the case of a 
first-order transition or continuous for critical wetting. 

From a theoretical point of view, it is much more con- 
venient to take the adsorbed film thickness, £, rather than 
the contact angle, as an order parameter for wetting tran- 
sitions and related phenomena. An interfacial Hamilto- 
nian is then minimized with respect to £ as is typically 
the case with the (mesoscopic) Landau-type field theo- 
ries and (microscopic) density functional theory (DFT) 
- where i can be easily determined from the Gibbs ad- 
sorption, a direct output of DFT. 

In this study, we examine the wetting properties of a 
simple fluid in contact with a spherical attractive wall 
by using an intermolecular interaction model with fluid- 
fluid and fluid-wall long-range forces. The curved geom- 
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etry of the system prohibits a macroscopic growth of the 
adsorbed layer (and thus complete wetting), since the 
free energy contribution due to the liquid-gas interface 
increases with the film thickness £, and thus for a given 
radius of a spherical substrate there must be a maximum 
finite value of ^ [ESQ- For the mesoscopic approaches, 
the radius of the wall i?, is a new field variable that in- 
troduces one additional ^-dependent term to the effec- 
tive interface Hamiltonian of the system, compared to 
the planar geometry, where the only ^-dependent term is 
the binding potential between the wall-liquid and liquid- 
gas interfaces. Furthermore, for a fluid model exhibiting 
a gas-liquid phase transition, such as ours, it has been 
found that two regimes of the interfacial behavior should 
be distinguished: R > Rc, in which case the surface 
tension can be expanded in integer powers of R~^ and 
R < Rc, where the interfacial quantities exhibit a non- 
analytic behavior Q. Moreover, for an intermolecular in- 
teraction model with fiuid-fluid long-range interactions, 
there is an additional R~^ log R contribution to the sur- 
face tension in the R > Rc regime These striking 
observations actually challenge all curvature expansion 
approaches. In addition, a certain equivalence between 
a system of a saturated fluid on a spherical wall and a 
system of an unsaturated fluid on a planar wall above 
the wetting temperature has been found [F, Some- 
what surprisingly, DFT computations confirmed this cor- 
respondence at the level of the density profiles down to 
unexpectedly small radii of the wall [g. 

Most of these conjectures follow from the so-called 
sharp-kink approximation (SKA) based on a simple 
piece- wise constant approximation of a one-body density 
distribution of the fluid, i.e. a coarse-grained approach 
providing a link between mesoscopic Hamiltonian the- 
ories and microscopic DFT. The simple mathematical 
form of SKA has motivated many theoretical investiga- 
tions of wetting phenomena as it makes them analytically 
tractable. At the same time SKA appears to capture 
much of the underlying fundamental physics for planar 
substrates (often in conjugation with exact statistical me- 
chanical sum rules [§]). 

However, as we show in this work, SKA is only qualita- 
tive for spherical substrates, even though the functional 
form of the fllm growth can still be successfully inferred 
from the theory Q. We attribute this to the particu- 
lar approximation of the liquid-gas interface adapted by 
SKA. In particular, since the ^-dependent contribution to 
the interface Hamiltonian due to the curvature is propor- 
tional to the liquid-gas surface tension, the latter plays 
an important role compared to the planar geometry. 

More speciflcally, the curved geometry induces a 
Laplace pressure whose value depends on both fllm thick- 
ness and the surface tension and so the two quantities 
are now coupled, in contrast with the planar geometry 
where a parallel shift of the liquid-gas dividing surface 
does not influence the surface contribution to the free 
energy of the system. We further employ an alternative 
coarse-grained approach, a modiflcation of the one origi- 



nally proposed by Napiorowski and Dietrich [10[ for the 
planar geometry, which replaces the jump in the density 
proflle at the liquid-gas interface of SKA by a continuous 
function restricted by several reasonable constraints. We 
show that in this "soft-interface approximation" (SIA) 
the leading curvature correction to the liquid-gas surface 
tension is 0(i?~^), rather than 0(i?~^ log R), in line with 
the Tolman theory. Once a particular approximation for 
the liquid-gas interface is taken, the corresponding Tol- 
man length can be easily determined. Apart from this, 
we flnd that the finite width of the liquid-gas interface 
significantly improves the prediction of the correspond- 
ing surface tension when compared with the microscopic 
DFT computations, which consequently markedly im- 
proves the estimation of the film thickness in a spherical 
geometry. 

In Sec. |lT]we describe our microscopic model and the 
corresponding DFT formalism. In Sec. IIIII we present 
results of wetting phenomena on a planar wall obtained 
from our DFT based on a continuation scheme that al- 
lows us to trace metastable and unstable solutions. The 
results are compared with the analytical prediction as 
given by a minimization of the interface Hamiltonian 
based on SKA. We also make a connection between the 
two approaches by introducing the microscopic model 
into the interfacial Hamiltonian. In Sec. IIVI we turn 
our attention to the main part of our study, a thin liq- 
uid film on a spherical wall. We show that the SKA 
does not perform as well as might be desired, in partic- 
ular, it does not account for a quantitative description 
of the liquid-gas surface tension which plays a significant 
role when the substrate geometry is curved. We then 
introduce SIA and present an asymptotic analysis with 
the new approach. Comparison with DFT computations 
reveals a substantial improvement of the resulting inter- 
face Hamiltonian even for very simple approximations of 
the density distribution at the liquid-vapour interface, 
indicating the significance of a non-zero width of the in- 
terface. We conclude in Sec. |V] with a summary of our 
results and discussion. Appendix A describes the contin- 
uation method we developed for the numerical solution 
of DFT. In Appendix B we show derivations of the sur- 
face tension and the binding potential for both a planar 
and a spherical geometry within SKA. Finally, Appendix 
C shows derivations of the above quantities, including 
Tolman's length, using SIA. 



II. DFT 

A. General formalism 

DFT is based on Mermin's proof [ll[ that the free en- 
ergy of an inhomogeneous system at equilibrium can be 
expressed as a functional of an ensemble averaged one- 
body density, p{r) (see e.g. Ref. for more details). 
Thus, the free-energy functional T[p] contains all the 
equilibrium physics of the system under consideration. 



3 



Clearly, for a 3D fluid model one has to resort to an 
approximative functional. Here we adopt a simple but 
rather well established local density approximation 



np] = J /hs (p(r))p(r) dr+ 



p(r)p(r')</.(|r-r'|)dr'dr, (1) 



where /Hs(p(r)) is the free energy per particle of the 
hard-sphere fluid (accurately described by the Carnahan- 
Starling equation of state), including the ideal gas contri- 
bution. The contribution due to the long-range van der 
Waals forces is included in the mean-field manner. To 
be specific, we consider a full Lennard- Jones 12-6 (LJ) 
potential to model the fiuid-fluid attraction according to 
the Barker-Henderson perturbative scheme 



0(r) = 





4e 



r < a 
r > a 



(2) 



where for the sake of simplicity the Lennard- Jones pa- 
rameter a is taken equal to the hard-sphere diameter. 

The free-energy functional, T[p\, describes the intrin- 
sic properties of a given fiuid. The total free energy in- 
cluding also a contribution of the external field is related 
to the grand potential functional through the Legendre 
transform 



n[p]=T[p]+ / pir){V{r)~p)dr, 



(3) 



where /i is the chemical potential and V (r) is the external 
field due to the presence of a wall W cM.^, 



Viv) 



PwJw4'w{\r~r'\) dr' elsewhere. 



(4) 



consisting of the atoms interacting with the fluid particles 
via the Lennard- Jones potential, (j)^) (r), with the param- 
eters aw and e^, and uniformly distributed throughout 
the wall with a density p^: 



(5) 



Applying the variational principle to the grand poten- 
tial functional, Eq. ([3]), we attain the Euler-Lagrange 
equation: 



Sp{r) 



p{r')(t){\r-r'\)dr' + V{r) 



0, (6) 



where J^hs[p] denotes the first term in the right-hand- 
side of (dl. In general, the solution of ^ comprises all 
extremes of the grand potential ^[p] as given by and 
not just the global minimum corresponding to the equi- 
librium state. Here we develop a pseudo arc-length con- 
tinuation scheme for the numerical computation of ^ 



that enables us to capture both locally stable and un- 
stable solutions and thus to construct the entire bifur- 
cation diagrams for the isotherms (details of the scheme 
are given in Appendix [X]) . 

The excess part of the grand potential functional ([3]) 
over the bulk may be expressed in the form 



n,,[p{r)]=~ J {p{p{r))~p{pb))dr+ 

+ 111 P W (P (r') - P (r)) <^(|r' - r|)dr'dr+ 
+ f p{r)V{r)dr, (7) 



where pb is the density of the bulk phase and 

~P {p) = p/hs (p) + ap^ - pp, (8) 

is the negative pressure, or grand potential per unit 
volume, of a system with uniform density p and a = 
\ J (j){\r\)dr = —^Trea^. In particular, the equilibrium 
value of the excess grand potential ([7]) per unit area of a 
two-phase system of liquid and vapour in the absence of 
an external field, yields the surface tension between the 
coexisting phases, 7;^. The prediction of jig as given by 
minimization of Q agrees fairly well with both compu- 
tations and experimental data as shown in Fig. [1] 




FIG. 1. Plots of surface tension as a function of dimension- 
less temperature, T/Tc. Solid line: numerical DFT results of 
our model scaled with e/fcs — 119.8K and a = SAA; trian- 
gles: computational results by Toxvaerd for a 12-6 LJ fluid 
using the Barker- Henderson perturbation theory [13 ] with the 
Percus-Yevick solution [l3| for the hard-sphere reference fluid 
and using the exact hard sphere diamet er JlSl] : circles: Monte 
Carlo simulations by Lee and Barker uM: squares: experi- 
mental results for Argon by Guggenheim [13]; dashed line: fit 
of experimental results to equation 7(r) = 70 (1 — T/TcY^^ 
by Guggenheim [T?]. The resulting coefficients are 70 = 
36.31dyn/cm and = f • 
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B. Translational symmetry: planar wall 

If the general formalism outlined above is applied on a 
particular external field attaining a certain symmetry, it 
will adopt a significantly simpler form. In the next sub- 
section we will formulate the basic equations resulting 
from the equilibrium conditions obtained from the mini- 
mization of ([71), for a spherical model of the external field, 
i.e. a system with rotational symmetry. But prior to 
that, it is instructive to discuss the zero-curvature limit 
of the above model, corresponding to an adsorbed LJ 
fluid on a planar wall, a system with translational sym- 
metry. 

For a planar substrate W — M.'^ x in Cartesian 
coordinates, the density profile is only a function of z, so 
that the Euler-Lagrange equation reads 

MHS(P(2))+ / p(z')*Pla(k-z'|)rfz'+ (9) 

+Voo{z)-fi = (VzeM+), 

where /ihs(p) = ^U^s(p)p) ^j^g chemical potential of 
the hard-sphere system. A fluid particle at a distance z 
from the wall experiences the wall potential: 



Voo{z)^Pu, MV^'^ + 2/" + (2 - z')^)dx'dy'dz' 
Jw 

z < 



oo 



-\^^P^e^ol[l,{^Y-\{^f) z>0. 

(10) 

^pia{z) in Eq. ^ is the surface potential exerted by the 
fluid particles uniformly distributed (with a unit density) 
over the x-y plane at distance z: 



(z) — II 4> {^^/ x'^ +y'^ + z^^ dydx 

(j) (^\/ + r^^ rdr (11) 



6 2/1 



5 ( 2.)^ _ 1 ( 

3 V z y 3 



f)- z>a. 



In the framework of DFT, the natural order parameter 
for the wetting transitions is the Gibbs adsorption per 
unit area: 

/>oo 

roo[p(^)] = / {p{z)-pb)dz. (12) 

Jo 



C. Rotational symmetry: spherical wall 

If the external field is induced by a spherical wall, W — 
{r e K'^ : r = |r| < i?}, the variational principle yields 

/•oo 

PHs{p{r))+ p{r )<^>sph {r,r')dr'+ (13) 
Jr 



where <&sph (^, ^') is the surface interaction potential per 
unit density generated by fluid particles uniformly dis- 
tributed on the surface of the sphere Br' centered at the 
origin at distance r, 

$sph(r,r')= / <l>i\r-r\)dr. (14) 
($Pia(|r-r'|)-$pia(|r + r'|)) 



r 
r 



(see also Appendix IB 1[) . The wall potential in Eq. 
for the spherical wall ly = {r e : |r| < R} is: 



3r 



4 
30 



9i? 



9R 



[r + i?)9 
r + 3R 



i?)9 



r-SR 



(r-i?)3 (r + i?)3 



(15) 



Replacing the distance from the origin r by the radial 
distance from the wall f = r — R, one can easily see that 
the external potential ([15]) reduces to the planar wall 
potential ([TU)) . for i? — >■ 00. Analogously to the planar 
case, we define the adsorption Tn as the excess number 
of particles of the system with respect to the surface of 
the wall: 



+VR{r) -p^O, (Vr > i?) , 



rfl[p(r)]= / (^) {p{r)-p,)dr. (16) 



III. WETTING ON A PLANAR SUBSTRATE 

In this section we make a comparison between the nu- 
merical solution of DFT and the prediction given by the 
effective interfacial Hamiltonian according to SKA for 
the first-order wetting transition on the planar substrate. 
We consider a planar semi-infinite wall interacting with 
the fiuid according to ([T0|) with the typical parameters 
Pw^w = 0.8e/a^ and a^i = 1.25(t that correspond to 
the class of intermediate-substrate systems [31 for which 
prewetting phase transitions can be observed. We note 
that wetting on planar and spherical walls is a multipara- 
metric problem and hence a full parametric study of the 
global phase diagram is a difficult task, beyond the scope 
of this paper. 



A. Numerical DFT results of wetting on a planar 
wall 

Figure [2] depicts the surface-phase diagram of the con- 
sidered model in the (A/i,T) plane, where Afi = p — psat 
is the departure of the chemical potential from its sat- 
uration value. The first-order wetting transition takes 
place at wetting temperature ksTw = 0.621e, well bellow 
the critical temperature of the bulk fluid UbTc = 1.006e 
for our model. The prewetting line connects the satura- 
tion line at the wetting temperature and terminates 
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at the prewetting critical point, ksTp^c = 0.724£. The 
slope of the prewetting line is governed by a Clapeyron- 
type equation [l^, which, in particular, states that the 
prewetting line approaches the saturation line tangen- 
tially at with 



dT 



= 0, 



(17) 



in line with our numerical computations. Schick and Ta- 
borek 20] later showed that the prewetting line scales as 
— A/i ^ (T — Tyj)'^^^. In Ref. [21], this power law was 
confirmed experimentally, such that 



knT,,, 



C 



Tn, 



3/2 



(18) 



with C w i 



211. A fit of our DFT resuhs with ^ leads 
to a coefficient C = 0.77, in a reasonable agreement with 
the experimental data - see Fig. [21 

Figure |3] depicts the adsorption isotherm in terms of 
the thickness of the adsorbed liquid film £ as a function of 
A/i for the temperature ksT — 0.7e and in the interval 
between the wetting temperature T^, and the prewetting 
critical temperature Tp^c- ^ can be associated with the 
Gibbs adsorption through 



Ap 



(19) 



for both finite and infinite R, where Ap — Pi"'* — p^"* 
is the difference between the liquid and gas densities at 
saturation. 

The isotherm exhibits a van der Waals loop with two 
turning points depicted as B and C demarcating the un- 
stable branch. Points A and D indicate the equilibrium 
between thin and thick layers, corresponding to a point 
on the prewetting line in Fig. [2] The location of the 
equilibrium points can be obtained from a Maxwell con- 
struction. Details of the numerical scheme we developed 
for tracing the adsorption isotherms are given in Ap- 
pendix |X1 



B. SKA for a planar wall 

For the sake of clarity and completeness we briefly re- 
view the main features of SKA for a planar geometry 
(details are given in Ref. [3]). 

Let us consider a liquid film of a thickness £ adsorbed 
on a planar wall. According to the SKA the density dis- 
tribution is approximated by a piecewise constant func- 
tion 



z<5, 

pI s < z < 

Pa Z > £, 



(20) 
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FIG. 2. The upper plot shows the deviation of the chemi- 
cal potential from its saturation value at prewetting (crosses) , 
and at the left (open squares) and right (filled squares) sad- 
dle nodes of bifurcation as a function of temperature. The 
dashed line marks the locus of the chemical potential at sat- 
uration for the given temperature, Ap — 0. The solid line 
is a fit to -App^{T)/{kBT^) = C((r - T^)/T^f/^ where 
the wetting temperature is ksT^ ~ 0.62l£ and the prewet- 
ting critical temperature is ftsTpmc = 0.724e. The resulting 
coefficient is C = 0.77. The lower plot shows scaled prewet- 
ting phase diagrams for different systems. The circles are 
DFT calculations for an attractive wall with CTu, — 1.25a and 
pw£w ~ 0.8e/cr^ (open circles) and pwEw ~ 0.75e/(T^ (filled 
circles). Experimental data ]2^: filled squares, methanol on 
cyclohexane (Kellay et. al. 1993) ,2^1]; open triangles, H2 
on rubidium (Mistura et al. 1994) 23]; filled triangles, He 
on caesium (Rutledge et al, 1997) i2J]; open squares, H2 on 
Caesium (Ross et al, 1997) [H]. 



where pg is the density of the gas reservoir and pj' is the 
density of the metastable liquid at the same thermody- 
namic conditions stabilized by the presence of the planar 
wall, Eq. PH)) and S ~ ^{a + aw). The off-coexistence of 



the two phases induces the pressure difference 
P^(Ai) -p(m) ~ ApA/i, 



(21) 



where p+ is the pressure of the metastable liquid and p 
is the pressure of the gas reservoir, and where we assume 
that Ap = p — Psat < is small. 
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FIG. 3. The upper graph depicts an ^-A/i bifurcation dia- 
gram for fcsT = 0.7e for a wall with pwEw = O.Se/o""^ and 
= 1.25(7. A/i is the deviation of the chemical poten- 
tial from its saturation value, ^sat- The prewetting transi- 
tion, marked by the dashed line, occurs at chemical potential 
^f^pw ~ — 0.022e. The inset subplots show the density pa'^ 
as a function of the distance z/a from the wall. The lower 
graph shows the excess grand potential ilex/£ as a function 
of A/i/e in the vicinity of the prewetting transition. 



The excess grand potential per unit area A of the sys- 
tem then can be expressed in terms of macroscopic quan- 
tities as a function of £ 



A 



(22) 



where 7^^^ and 7;^^"^ are the SKA to the wall-liquid and 
the liquid-gas surface tensions, respectively, and 'w^^^{£) 
is the effective potential between the two interfaces (bind- 
ing potential). In the following, we will suppress the ex- 
plicit /i-dependence. 



The link with the microscopic theory can be made, if 
the contributions in the right-hand-side of Eq. (1^^ are 
expressed in terms of our molecular model, which, when 
summed up, give the excess grand potential ([7]) where we 
have substituted the ansatz (1201): 



Iwl 



+2 .0 .00 

-Hr / / ^PiA\^ -^'\)dz'dz+ (23) 

rOO 

+ pt Voo{z)dz 



SKA 
llg 9 



2 /-O roo 



<^m!.{\z- z'\)dz'dz 



4 



(24) 



„SKA 



{£) ^Ap(p+ / <S>pUz')dz'dz- 

\ J 1-5 J z 

V^{z)dz\ (25) 



127r£2 



where we considered the distinguished limit 5 <^ i. A\s 
the Hamaker constant corresponding to the net interac- 
tion of an atom and the semi-infinite wall: 



(26) 



We note that the Hamaker constant is implicitly temper- 
ature dependent and that the attractive contribution of 
the potential of the wall enables the Hamaker constant to 
change its sign. Hence, in contrast with the adsorption 
on a hard wall, where the Hamaker constant is always 
negative, there may be a temperature below which its 
sign is positive (large pi) and negative above. Clearly, 
complete wetting is only possible for A < 0. 

Making use of only the leading-order term in ((25|) the 
minimization of (|22[) with respect to £ gives: 



ApA^ 



A 



0. 



(27) 



Hence, at this level of approximation the equilibrium 
thickness of the liquid film is: 



£eq 

When substituted into 
to leading order reads: 



/ A 



(28) 



\ 1/3 

l^67rA/9A/iy 

the wall-gas surface tension 



^SKA 
Iwg 



7,r 



9A 



1/3 



lApA^il 



2/3 



(29) 



7 



Equation (1281) can be confirmed by a comparison 
against the numerical DFT, see Fig. S] We observe that 
the prediction of SKA becomes rehable for | A/i| < O.Ole 
corresponding to a somewhat surprisingly small value of 
the liquid film, i « 5(7. Beyond this value, the coarse- 
grained approach looses its validity and also the prewet- 
ting transition is approached, both of which cause the 
curve in Fig. |3]to bend (see also Fig. [3]). It is worth 
noting that the only term in (p2)) having an ^-dependence 
and thus governing the wetting behavior, is the term re- 
lated to the undersaturation pressure and the binding 
potential, w^^^{C.). Clearly, 7;^ does not come into play 
in the planar case since the translation of the liquid-gas 
interface along the z axis does not change the free energy 
of the system. The situation becomes qualitatively differ- 
ent if the substrate is curved. Nevertheless, at this stage 
we conclude in line with the earlier studies, that SKA 
provides a fully satisfactory approach to the first-order 
wetting transition on a planar wall. 



b 




FIG. 4. Log-log plot of the film thickness as a function of 
deviation of the chemical potential from saturation, A/^, for 
ksT — 0.7e and wall parameters pm£w = O.Se/cr'^, a^y ~ 
1.25(T. The crosses are results from DFT computations. The 
solid line is the analytical prediction in Eq. (|27[) obtained from 
SKA. 



The corresponding excess grand potential now reads 



-AfiAp 



'SR^ 



R 



(31) 



where R — R + d. Within this approximation, the liquid- 
vapour surface tension becomes (see also Appendix [B|) 



(i?) 



(ex.) 



1 - 



2 \n{R/cr) 
9 {R/af 



O 



[{a/RY 



(32) 



and an analogous expansion holds for 7f,f^(-R). The 
^{R/iy^ correction to 7;g^^(oo) is due to the decay 
of our model. We note that short-range potentials lead 
to different curvature dependence of the surface tension, 
a point that has been discussed in detail in Refs. 0, H, 
l26l |. Interestingly, the 0{a/R) correction to the surface 
tension, as one would expect from the Tolman theory [13] , 
is missing. It corresponds to a vanishing Tolman length 
within SKA, as we will explicitly show in the following 
section. Although the value of the Tolman length is still 
a subject of some controversy, it is most likely that its 
value is non-zero, unless the system is symmetric under 
interchange between the two coexisting phases [1^ . This 
observation has been confirmed numerically in Ref. Q 
from a fit of DFT results for the wall-gas surface tension 
in a non-drying regime for the hard-wall substrate. Thus, 
the linear term was included by hand into the expansion 

§21 i- 

Finally, the binding potential within the SKA for the 
spherical wall yields 



(£;i?) 



„SKA , 



?;oo) 1 + 



R 



(33) 



where terms O {{6 / if , 5 / R.^^f^^ 
glected. 



have been ne- 



IV. WETTING ON A CURVED SUBSTRATE 



B. SI A for the spherical wall 



A. SKA for the spherical wall 

For the spherical geometry, SKA adopts the following 
form: 

r < i? + (5, 
p^^tir) = <( pt R + S <r < R + £, (30) 
Pg R + i < r < 00 . 



As an alternative to SKA, Napiorkowski and Diet- 
rich [13] proposed a modified version of the effective 
Hamiltonian, in which the liquid-gas interface was ap- 
proximated in a less crude way by a continuous mono- 
tonic function, the SIA. Applied for the second-order wet- 
ting transition on a planar wall, SIA merely confirmed 
that SKA provides a reliable prediction for such a sys- 
tem, in complete agreement with SIA. Formulated now 
for the spherical case, the density profile of the fiuid takes 
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FIG. 5. Sketch of the density profile according to SIA for a 
certain film thickness I. A piecewise function approximation 
is employed so that except for the interval (_R + £ — x/2, R + 
I + x/2) the density is assumed to be piecewise constant. 



the form: 



(34) 



= < 





Pt 

Pig{r- 
Pg 



R-l) 



r < R + 6 
R + 6 <r <R 
R + i 
R + i 



e- 



- I < r < i? - 

+ ^ < 7' < OO 



X 



Thus, a non-zero width of the liquid- vapour interface, Xi 
is introduced as an additional parameter. The density 
profile pig{-) in this region is not specified, but the fol- 
lowing constraints are imposed: 



Pig 



(-1) 



and 



Pig 



(f) 



Pa^ 



(35) 



with an additional assumption of a monotonic behaviour 
of the function pig{r). An illustrative example of pj}_j{r) 
is given in Fig. [5] The corresponding excess grand po- 
tential takes the form 



= - A^Ap 



[R + if-R^ , 



3R^ 



^V^r) 7/r(^ + ^) + ^ (^'^)' (36) 

taking i? -f ^ as the Gibbs dividing surface (so that £ 
is a measure of the number of particles adsorbed at the 
wall). The binding potential (see also Appendix IC 3[) is 
obtained from 



w""'^ {R,£) 



^-Pi I (P; - W (^) dr- 

2 



{pt-pfi{r))Vn{r)[^ dr, (37) 



where 4'fl(r) — ^Sphir, r')dr' - see Appendix IB II for 
the explicit form of the last expression. 

The wall-liquid surface tension remains unchanged 
compared to that obtained from SKA, Eq. (|24l) . How- 
ever, the liquid-gas surface tension now reads (see Ap- 
pendix [CT|) 



{R) = - 
1 

+ - 



«+x/2 ,^.2 

iPiPigir)) - p{Pref{z))) [ — ) dr- 

i?-x/2 ^-K^ 



Plg{r) {Plg{r') - Plg{r)) X 



X «'sph(?',r') dr'dr. 



(38) 



From now on, we neglect the curvature dependence 
of X s-nd Pig (•), as they would introduce higher-order 
corrections not affecting the asymptotic results at our 
level of approximation. This is also in line with previous 
studies which show that the Tolman length only depends 
on the density profile in the planar limit [l^. Then ([38)1 
can be written as 



(R) 



(oo) 



1 - 



2S„ 



R 



( HR/a) 



(39) 



where pref{z) — pfQ[—z) + pg<d{z) and S^o is the Tol- 
man length of the liquid-gas surface tension, as given by 
(Appendix [C2|): 



x/2 



{p(pig{z)) - Pre}) Zdz. (40) 



7fg (oo) y_x/2 

Here, p^e/ is the pressure at saturation, so that the Tol- 
man length is independent of the choice of the dividing 
surface. An immediate consequence of Eq. (PD|) is that 
within SKA the Tolman length vanishes, since pref{z) 
corresponds to the SKA of a liquid-gas interface. 

The equilibrium film thickness then follows from set- 
ting the derivative of (|36|) w.r.t. t equal to zero: 



1 dr^ei 
47ri?2 dt 



-AfiAp 1 



(41) 



2 ^^SIA 



d£ 



rB.+e+x/2 2 

Pi I p[gir^R~£)^R+s{r)(-) dr- 

lR.+e-x/2 ^R^ 



p[^(r-R-£)Vnir)[^y 



rR+e+x/2 
lR+e-x/2 

The last two terms of (I^T]) are of the form 
rx/2 



dr. 



x/2 



pl{r)fiji{R + £ + r)dr, 



(42) 



with fiir) = p+^i,,+sir) (|) and ///(r) = Vnir) {j^y. 
Since pig{r) is monotonic, i.e. pj^ does not change sign. 
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the mean value theorem can be employed such that 

rx/2 

p'ig{r)fiji{R + l + r)clr = 



x/2 



(43) 



^Pfi.iiiR 



^iji), 



for some ^ui G (—x/2, x/2), where we made use of 
/ p'ig{r)dr = -Ap. Substituting gS]) into (gl]) and set- 
ting the resulting expression equal to zero, we obtain: 



Ap 



R + e 



pt^R+siR + i + Ci) 1 



Vr{R + £ + Cii) 1 




(44) 



So far, there is no approximation within SIA. Equa- 
tion (j44l) can be simplified by appropriately estimat- 
ing the values of the auxiliary parameters ^7 and ^//. 
To this end, we employ a simple linear approximation 
to the density profile at the liquid-gas interface, taking 
—p'ig{r)/Ap 1/x in (03]). Furthermore, we expand //,// 
in powers of £/R,<t/£ 



fi{R + 



2Tipfea^ 




fii [R + 



where we assumed the distinguished limits r,d,a <^ i <^ 
R. Inserting (gS]) and (gB)) into (gS]) yields for ^^r 

From (gU, we obtain to leading order. 



p+^R+siR + e + ^i) 1 



R + e 



(48) 



2^ + ef.^nf^ ^ f^'^ 



VR{R + e + iii) 1 



27r 



R + e 



(49) 



Finally, substituting (g5]) and (gS]) into (gj) we have to 
leading order: 



ApA^,-^^fl\^) 



A 



(50) 



and hence, to leading order the equilibrium wetting film 
thickness is: 



1/3 



nSIA 



A 



GnlApAp, 



27f;,i/i?) 



(51) 



We note that this asymptotic analysis can be extended 
beyond ([5T|). by including terms 0{S/e), 0{e/R) and 
0{{x/e)^). The latter occurs due to the "soft" treatment 
of the liquid- vapor interface and is thus not present in 
SKA. 

In Fig.[S]we compare two adsorption isotherms {kgT = 
0.7e) corresponding to wetting on a planar and a spher- 
ical wall {R = lOOcr). The two curves are mutually hor- 
izontally shifted by a practically constant value, in ac- 
cordance with Eq. (|50|) . This implies that the curve for 
the spherical wall crosses the saturation line — 
at a finite value of e, and eventually converges to the 
saturation line as A/i~^ from the right, thus the finite 
curvature prevents complete wetting. The horizontal 
shift corresponds to the Laplace pressure contribution, 
A/i = 27j^^^(oo)/ (ApR), as verified by comparison with 
the numerical DFT, Fig. [T] All these conclusions are in 
line with SKA. However, the difference between SKA and 
SIA consists in a different treatment of 7;g(oo), compare 
(|B4p and (|C2[) . This is quite obvious, since the softness 
of the interface infiuences the free energy required to in- 
crease the film thickness. We will discuss this point in 
more detail in the following section. 
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FIG. 6. Isotherms and density profiles for a planar wall 
(dashed lines) and a sphere with R — lOOcr (solid lines) at 
fcflT = 0.7e and with wall parameters, pw£w = 0.8e/cr^ and 
CT„ = 1.25ct. To directly compare the planar to the spheri- 
cal case, the film thickness instead of adsorption is used as a 
measure. The subplots in the inset depict the density pcr^ as 
a function of the distance from the wall z/a and {r — R)/a for 
the planar and the spherical cases, respectively. The points 
A and A' are at the prewetting transitions. Points B, B' and 
C, C' correspond to the same film thickness. B is at satura- 
tion whereas C is chosen such that the film thickness £ is 20a. 
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TABLE I. Planar surface tensions (|C2|) . Tolman lengths (|40|) 
and the corresponding parameters for temperature fcsT = 
0.7e according to a given auxiliary function approximating 
the density distribution of the vapour-liquid interface. The 
parameters are from auxiliary function minimization. The 
surface tension given by numerical DFT computations is 7;g = 
0.517e/(T^ and p = (p; + Pg)/2. Note that in the tanh-case, 
the interface width is implicitly determined by the steepness 
parameter a. 



Auxiliary function pig{z) 



7fq^(oo) argument 



P 
P 

I i 



^tanh {az/a) 



0.544e/cr^ X- 
0.532e/a2 X' 
0.524e/a2 a 



4.0cr -0.07(T 
5.4cr -0.09(T 
:0.66 -0.11a 



FIG. 7. Numerical verification of Eq. ((SOj. The film thick- 
ness £ is fixed and corresponds to the adsorption Fr = 
3.905/ct^. The solid line corresponds to the analytical result, 
Ap - 2jf}^{oo)/ (ApR) = Ce, where 7,5^^(00) = 0.524e/a'^, 
see Table]!] The symbols denote the numerical DFT results. 



C. Comparison of SKA and SI A 

We now examine the repercussions of the way the 
liquid-gas interface is treated on the prediction of wetting 
behaviour on a spherical surface. As already mentioned 
in Sec. lIVBl the linear correction in the curvature to the 
planar liquid-gas surface tension, ignored within SKA, is 
properly captured by SIA. Furthermore, the presence of 
the Laplace pressure suggests that the liquid-gas surface 
tension plays a strong part in the determination of the 
equilibrium film thickness. This contrasts to the case of 
a planar geometry, where the term associated with the 
liquid-gas surface tension has no impact on the equilib- 
rium configuration. 

To investigate this point in detail, we will first com- 
pare the approximations of jig as obtained by the two 
approaches. For this purpose, we start with SIA for a 
given parameterization of the liquid-gas interface. As 
shown in Table HI we employ linear, cubic and hyper- 
bolic tangent auxiliary functions, where the latter vio- 
lates condition ([55]) negligibly. The particular parame- 
ters are determined by minimization of a given function 
with respect to the corresponding parameters. In Table 
U we display the planar liquid-gas surface tension asso- 
ciated with a particular parameterization and the Tol- 
man length resulting from Eq. (|40p for the temperature 
ksT — 0.7e. In all three cases the surface tension is 
close to the one obtained from the numerical solution of 
DFT and also the predictions of the Tolman length are in 
a reasonable agreement with the most recent simulation 
results [29l - l3l| . with thermodynamic results as well 
as with results from the van der Waals square gradient 
theory 

It is reasonable to assume that from the set of the 




FIG. 8. Plot of a dimensionless planar liquid-gas surface 
tension for the liquid-gas interface approximation p{z) = 
'''^''^ — ^ taiih {az / a) for ksT = 0.7e as a function of the 
steepness parameter a. The upper dashed line is the surface 
tension obtained from SKA, whereas the lower dashed line 
displays the surface tension obtained from numerical DFT. 



considered auxiliary functions, the tanh-approximation 
is the most realistic one, although the numerical results 
as given in Table |T] suggest that it is mainly the finite 
width of the liquid-gas interface, rather than the approx- 
imation of the density profile at this region, that matters. 
To illustrate this, we show in Fig.|8]the dependence of the 
surface tension on the steepness parameter a, determin- 
ing the shape of the tanh function. Note that the limit 
a — > 00 corresponds to the surface tension as predicted by 
SKA, jf^^ = 1.060e/cr2, for fc^T = 0.7e. Such a value 
contrasts with the result of SIA, which corresponds to the 
minimum of the function, and yields 7;^^"^ = 0.524£/cr^, 
in much better agreement with the numerical solution of 
DFT,7g^J = 0.517e/a2. 

Asymptotic analysis of the film thickness in Eq. (|50p . 
reveals that the film thickness for large but finite R re- 
mains finite even at saturation with £ ^ R3 in line with 
earlier studies, e.g. Refs. [1,11]. From Eq. O one also 
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recognizes a strong dependence of £ on the planar liquid- 
gas surface tension. In Fig. ^ we present the SIA and 
SKA predictions of the dependence on £ as a function of 
the wall radius. The comparison with the numerical DFT 
results reveals that for large R SIA is clearly superior, re- 
flecting a more realistic estimation of the liquid-gas sur- 
face tension. For small values of R (and £) we observe 
a deviation between DFT and the SIA results. This in- 
dicates a limit of validity of our first-order analysis and 
the assumption of large film thicknesses. 



40 



b 

20 - 



10 - 



10 



.4 



R/a 



10' 



10" 



FIG. 9. Film thickness at saturation (Ajj, = 0) as a function 
of the wall radius. The symbols correspond to the numerical 
DFT results. The dashed line shows the prediction according 
to Eq. ini), where -yfg^ioo) = 0.524e/a'^ (see Table H]). The 
dash-dotted line corresponds to Eq. (|5ip where 'yf^^{oo) — 
1.060e/cr^ is used instead of 'yfg^{oQ). The wall parameters 
are pwSw ~ O.Se/a^ and a^j ~ 1.25a at fesT = 0.7e. 




FIG. 10. Density profiles of the fluid adsorbed at the spher- 
ical walls of radii R = 104. Icr (dashed) and R — 210. 6(t 
(dashed-dotted) in a saturated state and at the planar wall 
(solid line) in an undersaturated state, Afi = — 0.015e . The 
wall radii correspond to the equality 27^^ ao/-^ ~ "^pI^mI 



for j = SIA (dashed) and j = SKA (dashed-dotted). 
ksT = 0.7e. The wall parameters are pwSw = 0.8e/cr^ 
= 1.2o". 



For 
and 



with a view to verify the correctness of its functional 
form. Here, we show that the coarse-grained effective 
Hamiltonian approach is capable of a quantitatively reli- 
able prediction of the adsorption phenomena on a spheri- 
cal wall (for a sufficiently large R), if the restriction of the 
sharp liquid-gas interface is dropped. However, the price 
we have to pay for, is one more parameter (compared to 
SKA) that steps into the theory. 



The occurrence of the undersaturation pressure and 
the Laplace pressure on the left-hand-side of Eq. ([50)) 
suggests a certain equivalence between the two systems 
of a planar and a spherical symmetry once the sum of 
the two pressures is fixed. In Fig. [TU]we test this equiv- 
alence on the level of a density profile, where DFT re- 
sults corresponding to the planar and the spherical case 
are compared, such that Ap\Afi\ = 2jfg{oo)/R, with 
j — {SIA, SKA}. A high value of 7/3(00) as given by 
SKA must now be compensated by a fairly large R. As 
we have seen in Fig. [51 the high value of R means that 
the saturation line A/i = is crossed by the adsorp- 
tion isotherm at large £, in agreement with the result 
depicted in Fig. [Hi However, for a given R, I as ob- 
tained by SKA is underestimated, which follows from 
([5T|) with 7ig(oo) = 7;g^^(oo) which is also consistent 
with the physical observation that high surface tension 
inhibits growth of the liquid film. 

Note that these results are not in conflict with previous 
studies [Sj , where the SKA has been applied for drying on 
a spherical hard wall and very good agreement was ob- 
tained with DFT computations. This is because in Ref. 
[8] the "exact" (i.e. obtained from DFT computations) 
liquid-vapor surface tension was implemented into SKA 



V. SUMMARY AND CONCLUSIONS 

We have re-examined the properties of a well known 
coarse-grained interfacial Hamiltonian approach, origi- 
nally proposed by Dietrich [1] for the study of wetting 
phenomena on a planar substrate and based on a sharp- 
kink approximation (SKA). SKA relies on approximating 
the density profile by a piecewise constant function and 
has proved to provide significant insight into interfacial 
phenomena as it is mathematically tractable and gives 
reliable results for a wide spectrum of problems. This 
theory is phenomcnological in its origin, but a link with 
a microscopic density functional theory (DFT) can be 
made, which allows to express all the necessary quantities 
in terms of fluid-fluid and fluid-substrate interaction pa- 
rameters. Comparison with numerical DFT reveals that 
SKA provides a fully satisfactory approach to the theory 
of complete wetting on a planar surface. 

One of the aims of this study was to demonstrate that 
for a spherical geometry the prediction quality of SKA 
regarding interfacial properties and wetting characteris- 
tics is limited. More specifically, we demonstrated that 
SKA satisfactorily determines the functional form of the 
asymptotic behaviour of the film thickness for large radii 
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of the substrate but leads to a significant quantitative 
disagreement in the prediction of the adsorbed fihxi thick- 
ness when compared against numerical DFT. The source 
of the deviation is the presence of the Laplace pressure 
that is not quantitatively captured within the framework 
of SKA. This contribution originates in the dependence 
of the free energy of the the liquid-gas interface on a po- 
sition of a dividing surface, a property that is absent in 
the planar case. 

We then showed that the properties of the effective in- 
terfacial Hamiltonian approach can be substantially im- 
proved if SKA is replaced by a soft-interface approxima- 
tion (SIA) , where the assumption of the sharp liquid-gas 
interface is replaced by a less restrictive approximation 
in which the interface is treated as a continuous function 
of the density distribution. We demonstrated that SIA 
allows for mathematical scrutiny as it is still analytically 
tractable, e.g. it provides the curvature expansion of the 
surface tensions (non-analytic in the wall curvature) with 
the leading-order term proportional to a / R. Moreover, it 
allows to express the corresponding coefficient, the Tol- 
man length, in a fairly simple manner and the values it 
predicts for the Tolman length are in a reasonable agree- 
ment with the latest simulation results. 

This is in contrast with SKA, where the linear term in 
the surface tension expansion is missing, i.e. the Tolman 
length vanishes. This observation is in a full agreement 
with the conclusion of Fisher and Wortis [l^ , since SKA 
treats the fluid in a symmetric way, and thus the Tolman 
length must disappear as for the Ising-like models. In 
other words, according to SKA, the surface tension of a 
large drop is equivalent to the one of a bubble, provided 
the density profiles of the two systems are perfectly an- 
tisymmetric in the planar limit. This is no more true 
for SIA, due to the asymmetry of the "local" contribu- 
tions to the surface tension, i.e. the first term on the 
right-hand-side of Eq. ((55)) . 

Furthermore, comparison with our numerical DFT re- 
vealed that the SIA results of the film thickness as a 
function of the wall radius are offer a drastic improve- 
ment to the ones obtained from SKA. This follows from 
the fact that the surface tension of the planar liquid-gas 
interface according to SKA is overestimated, which in 
turn underestimates the interface growth. 

It should be emphasized that all the theoretical ap- 
proaches we have considered in this work are of a mean- 
field character, i.e. do not properly take into account the 
interfacial fluctuations (capillary waves) at the liquid-gas 
interface. However, for our fluid model of a power-law in- 
teraction these fluctuations are not expected to play any 
significant role, since the upper critical dimension asso- 
ciated with the considered system is d* = 2 [s^]. Nev- 
ertheless, what one has to take into account in order to 
obtain the correct critical behavior, is the broadening of 
the interface at the critical region. Evidently, this feature 
is not provided by SKA. Consequently, within SKA the 
liquid-gas surface tension vanishes as f = 1 — ^^[5]. In 
contrast, the SIA provides the expected mean- field be- 
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FIG. 11. Plot of liquid-gas surface tension vs. t = 1 — 
T/Tc- The squares are the result of SIA, where a simple 
linear interpolant has been used to model the interface density 
profile. The surface tension has been obtained by minimizing 
the grand potential with respect to the interface width x- The 
solid Une is a fit to 7ig(oo)a^/e = Ct^^^, where the resulting 
coefficient is C = 3.4. The inset shows a plot of the interface 
width x/c" over t. The solid line is a fit to x = Cxt'" , where 
Cx = 2.0 and a = 0.57. 



havior "fig(po) ~ ^2 ^ as it is able to capture the interface 
broadening near the critical point (see Fig. [TT|). 

The SIA developed here, can be naturally extended 
by "softening" the wall-liquid interface in an analogous 
way as done for the liquid-vapor interface. However, 
such a modification would have presumably only negli- 
gible impact on the prediction of the thickness of the 
adsorbed liquid film, since the contribution to the ex- 
cess free energy from the wall-liquid surface tension has 
no ^-dependence and the change of the binding poten- 
tial is expected to be small. On the other hand, it 
may be interesting to find the influence of this refine- 
ment on quantities such as the density profile at con- 
tact with the wall. For this purpose a non-local DFT 
(e.g. Rosenfeld's fundamental measure theory) would be 
needed though [stlsllset. 

We also note that despite our restriction to a model of 
spherical symmetry, our conclusions should be relevant 
for general curved geometries and should capture some 
of the qualitative aspects of wetting on non-planar sub- 
strates. Of particular interest would be the extension of 
this study to spatially heterogeneous, chemical or topo- 
graphical substrates. Such substrates have a significant 
effect on the wetting characteristics of the solid-liquid 
pair (e.g. Refs. |37H4^ ). 
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Appendix A: Numerical metliods 



where Sij denotes the Kroenecker-Delta, which we have 
used in order to consider the grid size at the boundaries. 

This system of equations is solved using a modified 
Newton method, where each step Ap is rescaled with a 
parameter A such that = + XAp is bounded 

in (0,6/7r) in order to avoid the singularity of Eq. ([8|). 
Note that we have made use of the vector notation 
P'= {po, ■ ■ ■ , Pn)^ ■ In each Newton step n, the linear sys- 
tem of equations 



J • Ap = g(p" 



(A2) 



has to be solved, where the elements of the Jacobian 
matrix J are given by: 



dpj 



(A3) 



A2 



5i] ■ Mhs(p«) + — (2 - (^jo - ■ *Pla(kj - z,\). 



For our computations we employ dimensionless values. 
We use (J and e, as the characteristic length and energy 
scales, respectively. 



1. Density profile 

To obtain the equilibrium density profiles, the ex- 
tremal conditions Eqs. (jH]) and (IT5|) for the planar and 
the spherical case, respectively, must be solved numeri- 
cally. As both cases are of dimension one, the same nu- 
merical method can be applied and we restrict ourselves 
to presenting the numerical method for the planar wall, 

The domain M is restricted to an interval of interest 
[zo, zn] with boundary conditions p{z) — for z < zq and 
p{z) = pg for z > zpf. zq £ (0, 1) is typically chosen to be 
0.6. This can be done due to the repulsive character of 
the wall. The interval [zq, zj^] is then divided in a uniform 
mesh, Zi — zq + i ■ Az with i — 0, . . . N, where Az = 
{zn — zq)/N is the grid size. Subsequently, the integral 
in Eq. ([9]) is discretized using a trapezoidal rule inside 
the domain [zo,2;jv], whereas the analytical expression 



*pia (z) = / $pia(z')rfz' 



/ 9 

(45 (7) ~ 6 (2) 



if Z < CT 
if Z > CT 



is used for the integral outside that interval. Hence, we 
obtain a system of A'^-l-l nonlinear equations with {p^, i — 
0, . . . , as unknowns, namely: 

5i(Po, ■ • ■ ,PAf)— ms {Pi) + Vac{zi) - p+ (Al) 

+ Pff^Pla {zn - Zi) + 

Az ! 

+ 51 ~ '^JO ~ ^3^^ Pj^'Pla {\Zj - Zi\) =0, 



2. Adsorption isotiierms 

Solving Eq. (jAip will only give one density profile p 
for each chemical potential /i. However, in the case of 
a prewetting transition, there can be multiple solutions 
for the same chemical potential. From these solutions, 
only one is stable, whereas the other solutions are meta- 
or unstable (see also Sec. IIII A[) . In order to compute 
the full bifurcation diagram of the set of density profiles 
over the chemical potential, a pseudo arc-length contin- 
uation scheme is developed similar to the one employed 
by Salinger and Frink [43|. 

More specifically, we introduce an arc-length 
parametrization such that (p{s),p{s)) with s G R 
is a connected set of solutions of condition (jAip and 
where we have included the chemical potential n as an 
additional variable: 



g(/i,p) =0. 



(A4) 



The main idea of the continuation scheme is to trace the 
set of solutions along the curve parametrized by s. 

Assume that a point (^",p") at position s" on the 
curve is given, where n is the step of the continua- 
tion scheme being solved for. First, the tangent vector 

^) position s" is computed. This is done by dif- 
ferentiating g(s):=g (/z(s), p(s)) with respect to s. From 
(|A4p . it is known that g is a constant equal to zero on 
the curve of solutions {p{s), p(s)). Hence, the differential 
^ vanishes: 



dg 

ds 



djJ, 

ds 
dp 
ds 



= 0, 



where J is the Jacobian as defined in (IA3I) and 



dgi 
dp. 



= -1 



dpg 

dp 



Pla 



{zN - Zi) . 



(A5) 



(A6) 
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such that, 




FIG. 12. Sketch of one iteration step of the continuation 
scheme, x" and x""^^ are consecutive points of the iteration, 
where, x — (p.,p). xt is the tangent vector in x". By fol- 
lowing the curve of solutions in the direction of the tangent 
vector, the pseudo arc-length continuation scheme is able to 
trace the curve of solutions through turning points with re- 
spect to the parameter /i. 



The second term takes into account that pg for the den- 
sity at z > z^v depends on the chemical potential. In 
our computations, we have approximated ^ by —1. 
(jASp is the defining equation for the tangent vector 

We remark that this homogeneous system of linear 
equations leaves one degree of freedom, as we only have 
N + \ equations, but iV -t- 2 variables, {^t,Pt)- An ad- 
ditional equation is then used to maintain the direction 
of the tangent vector on the curve of solutions: 



/ n-l f „n-l\ 



Pt 



where 1^ ^ (^pj i)"^ ^ is the tangent vector of the pre- 
vious iteration. 

In a second step, an additional equation for a point at 
the stepsize 6 away from (/i", p") and in the direction of 

the tangent vector ^ /i^^^ ^ is set up. For this 

purpose we introduce a scalar product, which takes into 
account the discretization of the density profile into N 
intervals of length Az: 



((Mi,Pi) I (M2,P2)):=AiiA*2 + 



(A7) 



A, 



N 



3=0 

The norm with respect to this scalar product is defined 
as: 



||(m,p)||:=((m,P)I(m,P))'/'. 



(A8) 



The curve of solutions (/i(s),p(s)) is now parameterized 
by the arc-length with respect to the norm given above, 



dp, dp 
ds ' ds 



ds = 9. 



(A9) 



Linearizing the norm around s" and making use of the 
approximate tangent vector [pttPt) at s", one obtains 

((m?, p^) I + e)- Pis-), Pis'- + e)- = 9, 

(AlO) 

where we have made use of the normalized tangent vector 
such that 



\\{pt.,Pt)\\ = 1. 



(All) 



Insert ing (^"+i,p"+i) for {p{s"^ + 0) , p{s"- + 9)) into 
(jA10|) leads to the additional equation for the next point 
on the curve of solutions: 



-\p"+i): = 



Kn {p" 

((m?,p?)I(m"+'-m",p"+' 



p")) 



(A12) 
0, 



For a geometric interpretation of Eq. (jA12[) see Fig. [T^l 
To obtain p"+i) , (Pl2)) is solved together with 

(|A4[) . This is done using a Newton scheme. In each 
Newton step, the following system of linear equations is 
solved: 



(P?)^ 



^ J 



( Ap- 
^ Ap" 



_ f Kn [p 



n.m ^n,m 

7 P 



(A13) 



where we are considering the n-th step of the continua- 
tion scheme and the m-th step of the Newton method, 



such that Ap 



m . ,,n,m-t-l 



p" 



and Ap" 



Furthermore, we have made use of 



Az 



Pt, 



(2 ~ 5jQ~ SjN) Pt,j- 



Finally, (|A13p is solved using a conjugate gradient 
method, where the Jacobian (|A3[) of the system is ap- 
proximated by introducing a cutoff of 5 molecular diam- 
eters for the intermolecular potential 'I'pia. 



Appendix B: Surface tension and binding potential 
in the sharp-kink approximation (SKA) 

1. Surface tension 

According to Gibbsian thermodynamics, the surface 
tension is the free energy cost to increase an interface 
by unit area, i.e. the excess free energy (excess grand 
potential for an open system) per unit area with respect 
to the corresponding uniform phases. Within SKA, the 
liquid- vapour surface tension can be obtained from ([7]), 
with 



P(r) 



PA, 
PB 



(Bl) 
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where n Vb = and V^i U Vb = R^. The convenience 
of the expression for the excess grand potential as given 
by ([7]) becomes evident now, as pA — Pi, Pb — Pg 
and no external field, only the second term in ([7|) mat- 
ters. One then gets an immediate result for the liquid-gas 
surface tension, 



7^^ 



where 



■^^^-^/(Va,Vb), (B2) 



; / / (^(|ri-r2|)dridr2. (B3) 



For the surface tension of a planar interface we have 

= V^<o and V_b = ^z>o such that 



$Sph(r-, r' 

r27T pTT 



I /.(r+r')^ 



= 27r— 
r 



/ / (t){\Y~-v'\)r'^sixid'dd'dif' 
Jo Jo 

4) (Vt) dt 



(r — r')2 
oo 



(r — r')^ 



- / (j}[y/i]dt 

(r+r'Y 



= -(<i>Pia(k-r'|)-<i>pia(|r + r'|)), 
and for r > R 



(B6) 



/(V.<o,V.>o) ^ 1 
A 2 



oo Jo 



*Pla(k-z'|)dz'dz, 



with 'I'pia defined by pTI) . Thus, for the liquid-gas sur- 
face tension we obtain: 



7i3(oo) 



Ap2 



2 ^0 /-oo 



oo Jo 



^P\i,{\z- z'\)dz'dz 



---nAp^ea^. 
4 



(B4) 



^R{r)= / ^sphl^J^Oc^^' 
Jo 



30 



r+9B. _ r-9R. 

r-3i?. _ _r+3R 



-a 



(B7) 



R^a <r 



2621 9_ 

15 cr 5ct2 



27 I r+9fl ^2_r+3_R ^ I? i ^ 



Note that expression (jB5p gives a vanishing Tolman's 
length. 



2. Binding potential 



In the case of a spherical symmetry, i.e. a drop of 
liquid of radius i?, Va = {r G : |r| < i?} and Vs = 
{r e M'' : |r| > i?} the surface tension becomes 



7zs(i?)=-Ap' 



47ri?2 



'7? Jo 

~Jr 



-^j <f>Sph (r, r') dr'dr 



=7;g(oo) ^1 



2 ln(i?/(7) 
9 (i?/a)2 



O ((a/i?)^ 



(B5) 



where Ap = pi-Pg and $sph (r, r') = /^^^^ (|r - r'|) dr' 
can be advantageously expressed in terms of $pia: 



The binding potential of a system possessing two in- 
terfaces is the surface free energy per unit area of the 
system minus the contribution due to the surface ten- 
sions of the two interfaces. It expresses an effective in- 
teraction between the interfaces induced by the attrac- 
tive forces. If, analogously to the analysis above, we de- 
fine three disjoint subspaces Yw, Va, and Vb, such that 
Yw U Va U Vb = the density distribution of the 
wall-liquid-gas system within SKA is 



0, reYw 
p(r) = { pi, vgYa , 
Pg r e Vb, 



(B8) 



which when substituted into (O gives for the excess grand 
potential: 

^-ApApVA - pi^I{Yw,Ya) - p1i{Yw,Yb) - 
-(Ap)2/(Va,¥b)+ / V^(r)p(r)dr. (B9) 
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We now rearrange the terms in (jB9p . such that 



is obtained by substituting (|C1|) into ([7]) with V{r) =0 



_ ^,,/s,„^A , „SKA I SKA , „.SKA 



(BIO) 

where A = ^q^w '^^ surface of the wall and A' = 

/a(v»-uVA) surface of the liquid-gas interface. 

We obtain 



1 



ill -Pi'/(Vw,VA U Vb) +Pi / l^(r)dr 



7f^ = -^(Ap)'/(VM^UV^,VB), 



(Bll) 
(B12) 



and the binding potential w^^^ involving the remaining 
contribution 



SKA 



(f) =^ (^piApI(Nw,^B) - ApJ^ V{r)dr 



w 



(B13) 

Having obtained the expressions of I{X, Y) for systems 
possessing translational or spherical symmetry, we can 
evaluate the binding potential in the planar case by mak- 
ing use of = R2 X (-00,(5], Va = x {6,e) and 
Vs =K^ X [£,oo): 



/ poo poo 




In the spherical case we make use of Yw ~ {r E M.^ : 
\r\ < R. + S}, Ya = {r e : R + S < \r\ < R + 1} and 
Vs = {r e . |j.| > ^ ^_ £} obtain 



e 



^SKA j^ypn^re ^ ^SKA ^) ^ _ ^ (3^5) 

where we have neglected terms O ,5/R, j^^^jr^ ■ 



(C2) 



A 

x/2 

{p{pig,oo{z)) - p{pref{z))) dz + 

-x/2 
00 pOO 

I Pig, 00 (z) {Plg,ao{z') Pig, 00 

(z)) X 

-00 J —00 



X ^p\a.iz, z')dz'dz, 



where pref{z) denotes the density of a given bulk phase, 
i.e. pref{z) = piO{—z) + pgQ{z) such that at saturation 
P{Pref{z)) = Pref = const. Wc notc that in the above 
approximation the contribution due to the excess local 
pressure is generally non-zero (in contrast to the SKA). 

In the spherical case, the density profile is 



Pi r<R-x/2 
PigAr)^{ Pia{r~R) |r - i?| < x/2 , (C3) 
Pg r>R + x/2 



and the surface tension of a liquid drop of radius R is 



47ri?2 



(C4) 



R-x/2 
00 pOO 



(piPigMr)) ~ Pref) (-^) dr+ 



"'0 



PlgA^) (PlgM^') - Plg-nir)) X 



X <&Sph(r,r') dr'dr. 



Appendix C: Surface tension, binding potential, and 
the Tolman length in the soft-interface 
approximation (SIA) 

1. Surface tension 

The surface tension of a planar liquid-gas interface in 
the soft-interface approximation 



PI z < -x/2 

Plg.,oc{z) = ^ pig{z) \z\ < x/2 
Pg Z> x/2 



(CI) 



2. Tolman length 



Here we calculate the Tolman length as given by SIA by 
a direct comparison of (|C2|) and (|C4|) . We first compare 
the second terms of (IC2p and (jC4p . For this purpose we 
define 



hnir, r') = p/g,_R(r) {pig,R{r') - p;g,fl(r)) 
and h{r, r') = pig,oo{r) (p;g,oo(r') - pig,oa{r)) (C5) 
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and taking use of (jB6[) we can express the double integral 
in (EH) as 



2 roo pOD 



hR{r,r')^Sphir,r') (^) dr'dr 
£ Jo Jo 

- / h{r,r'){^pU\r-r'\y 
£ J-rJ-r 

-$Pia(|2i? + r - r'D) (^1 + ^) (l + ^) dr'dr 

9 /•OO /"CXD 



£ J-Rj-R 



h{r, r')$pia(|r--r'|)x 



r \ 



a 

e 



2 /'OO /-OO 



/i(r, /)$pia(|r--r'|)x 



RJ-R. 



a 

e 



2 /"OO /'OO 



h{r, /)$pia(|r--r'|)x 



OO J — OO 



X 1 + 



Comparison with the double integral in ()C2p then yields 
— h{r, r')$Pia(k - r'\) (r + r') dr'dr+ 



O 



f\n{R/(7] 



a 

7r 



2 roo /'OO 

r (/i(r, r') + h{r' , r)) $Pia(|r - r'\)dr'dr+ 



— OO — OO 



7r 



2 roo poo 



r {pig,oo{r') - Pig,oo{r)) x 



OO <J — OO 



x<i>P,a(k-r'l)rf^'rfr- + 0(^i||^ 

In the following, we focus on the asymmetry of the model 
due to the contribution of the pressure but for simplic- 
ity we assume that the density profile is symmetric. In 
this case, the integrand in the above expression is anti- 
symmetric with respect to the reflection transformation 
r — >■ —r and r' — > — r' and the term 0{(j / R) vanishes. 

For the difference of the first terms of ((64| and ([C2|) 
we obtain 



2 poo 

{p{Plg,oo{z)) -Pref)dz 

— OO 

2 /-x/a 

{p{Plg.oo{z)) ~Pref)zdz + 

-x/2 

+ 0{ia/Rf) , 



2al 

eR 



(C7) 



yielding a Tolman length 
1 rx/z 



SIAr ^/ {p{pig{z))-pref)zdz. (C8) 

7;g looj J-x/2 



Note that in line with [28|, the Tolman length does not 
depend on the choice of the dividing surface. 



3. Binding potential 

The extension of the expression for the binding poten- 
tial, Eq. (jBlsp . as given by SKA is rather straightfor- 
ward. We consider the density distribution as follows 



P(r) = 








Pi 






r e Vab 


Pa 


r e Vb, 



(C9) 



for Nw a sphere of radius R + 5, Yw U Ya a sphere of 
radius R + £ — xf^, "^w U Ya U Yab a sphere of radius 
R + £ + x/2 and Vw' U Va U Vb U Yab = K^. Such a 
model is relevant for the study of wetting on a spherical 
(i? finite) and on a planar (i? — > oo) wall. It should be 
noted that in contrast to SKA, this density distribution 
is not piecewise constant, due to the position dependent 
part of p(r) in the region Yab- Furthermore, we define 
the following operators 

[XY] ^ -i ^ ^ (p (r) - p (r'))' (|r - r'|) dr'dr 
[XYU = -lf f {p^i{r)- p^,{r')f<j>{\v -r'\) dr'dr 

^ Jx JY 

[XY]i, ^-IJ^l {Pi9 (r) - Pig ir')f (|r - r'|) dr'dr , 



with (r) = piXR3\v„ir) and pig (r) = PiXYwUYAi^) + 
p/g(r)xv^B(r) -l-pgXVB(r), where xx{r) is the character- 
istic function of a subset X. Using this convention, the 
wall-liquid and liquid-gas surface tensions can be respec- 
tively expressed as 



l^oi =j ([VtvVa] + [Yw (Vab U Yb)U+ 



Pwiir)V{r)dr 

lig =j (^[VabVb] + ^[Yab^ab] + [V^ (V^B U Vb)] + 
+ [Vw (VabUVb)];^- 

(piPig i^)) - Pref) dr] , 
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where A = AttR^ . When this is subtracted from the In spherical coordinates, the binding potential reads: 
surface grand potential ([7]), which can be written as 



+ [Va(Vas U Vs)] + ]^[Yab'^ab] + [VabVb] 
{p{Pia {'c))-p{Pref (r)))rfr+ 



j p{v)V{v)dv^ , 



(CIO) 



one obtains for the binding potential: 



v^^^ = J {[^w (Vab U Vs)] - [Yw i^AB U Yb)U- 



[Vw (Vab U WB)]ig + V{r) (p (r) - p^i (r) dr) . 



R+5 i-oo 



w^'^= / / (^) pi{pi-pir'))x 

(r, r') dr'dr+ 

{p(r)~pi)VR{r)(^^ydr. (Cll) 
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